/ 


HjfU+y*  *  •> 

\  ‘  4 
A 


mm 


mb 

mat 

SKjS*; 
% 


file  COPJ 


AD- A 186  631 


2a.  StCUKIfY  CLASSIHLA I IUN  AU  I  HOKII  Y 


‘  DOCUMENTATION  PAGE 

I  lb  RESTRICTIVE  MARKINGS 


Form  Approved 
OMB  No.  0704-0188 


2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

ASE  87-312 


3  .  DISTRIBUTION  /AVAILABILITY  OF  REPORT 
Approved  for  public  release; 
distribution  unlimited 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

A*FOS*'T*'  87*1  468 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

Mississippi  State  University 


6b.  OFFICE  SYMBOL  I  7a.  NAME  OF  MONITORING  ORGANIZATION 


(If  applicable) 


6c.  ADDRESS  (City,  State,  and  ZIP  Code) 

Department  of  Aerospace  Engineering 
Mississippi  State,  MS  39762 


Air  Force  Office  of  Scientific  Research 


7b.  ADDRESS  (City,  State,  and  ZIP  Code) 

Directorate  of  Mathematical  &  Information 
Sciences,  Bolling  AFB,  D.C.  20332-6448 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 
AFOSR 


8b.  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable) 

NM  Grant  AF0SR-85-0143 


8c  ADDRESS  (City,  State,  and  ZIP  Code) 

Air  Force  Office  of  Scientific  Research 
Bolling  AFB,  D.C.  20332-6448 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM  PROJECT 

ELEMENT  NO.  NO. 

61102F  2304 


a/3 


I  WORK  UNIT 
lACCESSION  NO. 


COSATI  CODES 


SUB-GROUP 


18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 
?Grid  Generation,'  Curvilinear  Coordinates ,'  Numerical 
Methods j  Computational  Fluid  Dynamics  . 


19  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Numerical  solution  of  the  partial  differential  equations  for  the  generation  of  surface 
grids  requires  a  specification  of  the  forcing  function  which  depends  on  the  geometry 
of  the  given  surface.  The  data  for  the  surface  is  usually  in  discrete  form.  Methods 
have  been  developed  which  fit  a  function  over  the  given  data.  For  complicated  shapes, 
e.g.,  an  airplane,  the  functional  fit  and  the  eventual  grid  generation  for  the  fuselage 
and  wings  are  done  separately  and  then  integrated  later.  ( .  y  ,/■  ( $ 


.  OCT  2  6 1967 1 


20  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 
13  UNCLASSIFIED/UNLIMITED  □  SAME  AS  RPT 
22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 
Capt.  John  Thomas 

DO  Form  1473,  JUN  86 


□  OTIC  USERS 


21  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 

22b  TELEPHONE  (Include  Area  Code)  22c.  OFFICE  SYMBOL 

(202)  767-5026  NM  _ 


Previous  editions  are  obsolete 


87  10 


0  o  2 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 

UNCLASSIFIED 


1 1 .  TITLE  ( Include  Security  Classification) 

Generation  of  Surface  Grids  Through  Elliptic 

and  Missile  Configurations 

Partial  Differential  Equations  for  Aircraft 

12.  PERSONAL  AUTHOR(S) 

Z.U.A.  Wars! 

13a.  TYPE  OF  REPORT 

Interim 

13b.  TIME  COVERED 

from  4/86  to  6/87 

14.  DATE  OF  REPORT  (Year,  Month,  Day) 

1987  July  14 

15.  PAGE  COUNT 

25 

b: 


AFOSR-TK*  8  7-1  468 


Generation  of  Surface  Grids  Through  Elliptic  Partial 
Differential  Equations  for  Aircraft  and 
Missile  Configurations 


by 


Z.U.A.  Warsi 

Department  of  Aerospace  Engineering 
Mississippi  State  University 
Mississippi  State,  MS  39762 


Interim  Report 
April  1986  -  June  1987 
Submitted  to  the  Air  Force  Office  of  Scientific  Research 


|  Accession  For  j 

NTIS  GRA4I 
DTIC  TAB 
Unannounced 
Justifloatl* 

~ r 

□ 

an 

By 

Distributioi 

Availablllt 

1/ 

-y  Codes 

Olst 

M 

Avail 

Spec 

and/or 

lal 

Bolling  Air  Force  Base,  D.C.  20332 
Grant  No.  AFOSR-85-01 43 


1 .  INTRODUCTION 


The  problem  of  numerical  grid  generation  is  of  current  interest  in  many 
branches  of  engineering  particularly  in  aeronautics,  mechanical,  and  civil 
engineering.  The  spatial  grids  are  generated  either  by  algebraic  methods 
using  various  spline  and  transfinite  interpolations  or  by  solving  certain 
partial  differential  equations.  In  this  regard,  a  book  [ 1 ] ,  review  articles 
[2],  [3],  and  conference  proceedings  [4,  5,  6]  may  be  consulted. 

The  progress  report  under  consideration  is  in  reference  to  the  problem 
of  grid  generation  in  a  given  surface  by  solving  a  set  elliptic  partial 
differential  equations.  The  mathematical  model  used  in  the  present  research 
has  been  developed  by  the  author  under  previous  AFOSR  grants  and  has  been 
described  in  publications  [7  -  16].  Reference  [16]  has  been  attached  as  an 
appendix  in  which  Eqs.  (19)  and  (21)  describe  the  mathematical  model. 

The  developed  mathematical  model  has  been  used  to  generate  the  Car¬ 
tesian  coordinates  as  functions  of  the  curvilinear  coordinates  when  the 
surface  in  which  the  coordinates  are  to  be  generated  has  been  specified 
either  analytically  or  by  discrete  data  points.  In  most  practical  situa¬ 
tions  the  surface  is  usually  specified  by  discrete  data,  and  therefore, 
there  is  a  need  to  develop  computer  routines  to  fit  a  global  equation  of  the 
form  F(x,y,z)  -  0.  For  complicated  body  shapes  e.g.  an  airplane,  we  have  to 
divide  the  surface  into  suitable  sections  and  then  fit  the  functions  sepa¬ 
rately.  The  need  for  specifying  the  function  F  is  due  to  the  fact  that  the 
forcing  term  R  in  Eqs.  (19)  or  (21)  of  the  appendix  is  derived  from  F.  In 
the  period  under  consideration  we  have  devised  various  routines  for  the 
above  noted  purpose  and  they  are  described  in  Section  2.  The  developed 


computer  code  has  been  used  in  many  complicated  shaped  surfaces.  Figure  1 
shows  some  current  results. 

2.  NUMERICAL  SCHEMES 

A  computer  program  for  the  numerical  solution  of  Eq.  (21)  has  been 
developed  by  using  point  and  line  SOR.  The  essential  difference  between  the 
grid  generation  in  a  flat  space  and  in  a  curved  surface  is  in  the  appearance 
of  the  forcing  term  on  the  right  of  Eqs.  (19)  and  (21).  The  term  R  in  Eqs. 
(19)  and  (21)  depends  on  the  first  and  second  partial  derivatives  of  the 
fitted  function  F  with  respect  to  the  Cartesian  coordinates.  (Refer  to  [ 1 5 ] 
for  the  formula  of  R  in  terms  of  F).  The  following  routines  have  been 
developed  and  used  in  a  variety  of  problems. 

1.  Multidimensional  least  square  technique. 

2.  "Overlapping"  least  square  technique. 

3.  Fourier  decomposition  of  each  section  of  a  surface  and  then 
blending  them  to  obtain  the  equation  of  a  square. 

Beside  the  development  of  the  above  techniques,  we  have  also  developed 
the  following  two  separate  programs  in  the  solution  of  Eq.  (21). 

1.  For  the  acceleration  of  the  iterative  process,  a  "multigrid"  tech¬ 
nique  in  the  solution  of  Eq.  (21)  has  been  incorporated.  From  the 
test  cases  conducted  so  far  It  looks  that  the  solution  time  with 
multigrid  is  not  significantly  lower  than  that  of  SOR.  Research  in 
this  area  is,  however,  continuing. 

2.  A  routine  for  the  calculation  of  the  optimum  acceleration  parameter 
has  been  developed  which  works  quite  well  with  the  SOR  technique. 


3.  CONCLUSIONS 

Computer  programs  have  been  developed  which  solve  the  problem  of  nu¬ 
merical  coordinate  generation  in  a  given  surface.  The  mathematical  model  is 
formed  of  a  set  of  elliptic  equations  and  thus  the  resulting  coordinates  in 
the  surface  are  sufficiently  smooth.  The  grid  generation  in  arbitrary 
surface  in  accomplished  by  dividing  the  surface  into  suitable  subregions  and 
then  Joining  them  to  obtain  the  coordinates  in  the  integrated  surface. 
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ABSTRACT 

This  paper  is  devoted  to  an  analytical  comparison  of  the  various  elliptic  partial-dif- 
ferential-equation  (PDE)  models  which  are  in  current  use  for  grid  generation.  These 
comparisons,  particularly  between  the  equations  from  the  Laplace-Poisson  system  and 
the  equations  from  a  Gaussian  approach,  have  yielded  useful  expressions  connecting 
the  3D  Laplacians  and  the  surface  Bel  tram  ians.  This  effort  has  specifically  been 
successful  when  the  transverse  coordinate  leaving  the  surface  is  orthogonal  to  the 
surface.  Equations  which  are  derivable  from  Cartesian-type  Poisson  equations  and 
those  obtained  by  using  the  variational  principle  in  surface  coordinates  have  also  been 
considered. 


1.  INTRODUCTION 

The  problem  of  generating  spatial  coordinates  by  numerical  methods  is  of 
much  importance  in  many  branches  of  engineering  mechanics.  A  review  of 
various  methods  of  coordinate  generation  in  both  two-  and  three-dimensional 
Euclidean  spaces  is  available  in  Thompson,  Warsi,  and  Mastin  [1,  2]. 

This  paper  is  exclusively  directed  to  a  collection  and  analytical  compari¬ 
sons  of  the  various  elliptic  partial-differential-equation  (PDE)  models  which 
are  currently  in  use  for  numerical  coordinate  generation. 

The  theory  of  grid  generation  does  not  depend  on  any  set  of  so-called 
conservation  laws,  and  thus  a  variety  of  equations  and  methods  of  different 
characters  can  be  used  to  obtain  the  grids.  Any  cr  isistent  method,  depending 
either  on  the  solution  of  PDEs  or  any  algebraic  method,  can  be  used  to  obtain 
intersecting  trajectories  in  either  2D  or  3D  Euclidean  space.1 


•Research  supported  by  the  Air  Force  Office  of  Scientific  Research  under  Crant  AFOSR-85-0143. 
1  The  question  of  space  comes  into  the  picture  when  it  is  realized  that  all  the  memc  coefficients 
g,  cannot  be  selected  arbitrarily.  In  fact,  these  metric  coefficients  must  be  so  selected  that  a  set 
of  second-order  PDEs  become  satisfied.  For  details  refer  to  Warsi  (3). 
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In  light  of  the  above  discussion,  any  consistent  set  of  PDEs  is  sufficient  to 
form  a  mathematical  model.  In  bounded  domains  the  most  natural  choice  is 
;  that  of  a  mathematical  model  formed  of  elliptic  PDEs.  The  simplest  set  of 
equations,  which  also  yields  the  smoothest  grid,  consists  of  the  Laplace 
equations  of  the  curvilinear  coordinates  in  the  Cartesian  physical  space.  As  a 
next  logical  step,  a  set  of  Poisson  equations  can  be  selected  so  as  to  have  a 
degree  of  control  of  the  distribution  of  grid  lines.  On  inverting  these  equations 
a  set  of  quasilinear  PDEs  are  obtained  (see  Equations  (2),  (4)]. 

After  the  development  of  the  problem  of  grid  generation  through  solving 
the  inverted  forms  of  the  Laplace  and  Poisson  equations  in  2D  Euclidean 
domains  by  Allen  [4],  Winslow  [5],  Chu  [6],  and  Thompson  et  aL  [7],  a  logical 
extension  was  to  use  the  same  equations  for  3D  domains,  as  has  been  done  by 
Thompson  and  Mastin  [8].  Parallel  to  the  above-noted  developments,  Warsi 
[3,  9,  10]  proposed  a  Gaussian  approach  which  basically  generates  surfaces 
and  thus  can  be  used  either  for  generating  curvilinear  coordinates  in  a  given 
surface  [11]  or  for  generating  3D  coordinates  by  generating  a  series  of 
surfaces  starting  from  the  data  on  the  given  surfaces  [9,  12].  The  Gaussian 
approach  in  fact  depends  on  a  manipulation  of  the  formulae  of  Gauss  for  a 
surface,  and  thus  the  resulting  equations  have  the  surface  coordinates  as  the 
independent  variables.  This  manipulation  introduces  the  Beltramians  of  the 
curvilinear  coordinates  and  the  sum  of  the  principal  curvatures  of  the  surface 
in  a  natural  way.  Since  the  formulae  of  Gauss  for  a  surface  hold  true  for  any 
allowable  coordinate  system  introduced  in  the  surface,  the  equations  pro¬ 
posed  by  Warsi  [3,  9,  10]  must  also  have  the  same  properties.  Further, 
because  of  the  use  of  Gauss  formulae,  the  proposed  equations  are  optimal  in 
the  sense  that  the  number  of  terms  and  the  amount  of  information  in  the 
equations  is  just  sufficient  for  the  generation  of  either  surfaces  or  coordinates. 

In  this  regard  we  can  justifiably  call  a  mathematical  model  “optimal”  if  it 
can  be  reduced  to  the  form  of  the  proposed  equations  [cf.  Equation  (19a)]. 

Beside  the  Laplace-Poisson  system  and  the  Gaussian  approach,  we  have 
also  derived  the  surface  generating  equations  by  the  use  of  the  variational 
principle.  The  resulting  equations  are  near  optimal  for  the  Gaussian  system. 


2.  BASIC  ELLIPTIC  MODELS 

Poisson  Equations  as  Grid  Generators 

Since  the  publication  of  the  TTM  method  [7],  there  has  been  extensive  use 
of  the  Poisson  equations  in  the  physical  r-space  to  generate  both  2D  and  3D 
grids  [2].  In  practically  all  cases  the  main  aim  is  to  have  those  equations  in 
which  the  computational  coordinates  appear  as  the  independent  variables. 


Synopsis  of  Elliptic  PDE  Models  297 

and  therefore,  the  Poisson  equations  have  to  be  inverted  by  making  the 
physical  coordinates  r  =  (r,  y,  z)  the  dependent  variables.  This  inversion  can 
of  course  be  carried  out  in  a  nontensorial  manner  by  using  the  chain  rule  of 
partial  differentiation.  However,  it  is  much  simpler  to  follow  either  of  the 
following  two  methods  to  attain  the  same  result. 

Method  1:  Inner  multiplication  of  Equation  (A2)2  by  g';  and  the  use  of 
Equation  (A4)  results  in 


g"r.ji  +  (vV)rjk  =  0. 


(1) 


Introducing  the  second-order  differential  operator 


D 


=  ggri 


a2 

dx{dx> 


into  Equation  (1),  we  get 


Dr+  g(vzxk)r  k  =»  0,  (2) 


which  is  the  desired  equation  in  vector  form  in  the  computational  plane. 

Method  II:  In  Equation  (A3),  writing  4>  —  r  —  (x,y,  z)  and  again  using 
Equation  (A4),  we  obtain  the  same  equation  (2). 

To  form  a  closed  system  of  equations  from  Equation  (2)  one  has  to  specify 
the  Laplacians  V2xk  in  a  suitable  fashion.  As  described  in  Reference  [3],  a 
general  specification  for  the  Laplacians  is  to  write 


vV-g'fP*.  k  -1,2  or  k- 1,2,3,  (3) 


and  then  Equation  (1)  takes  the  form 

g^  +  P^rJ-O.  (4) 

In  Equation  (3),  the  P*  =  P*  are  intended  to  be  arbitrarily  specified  control 
functions.  Equating  the  right-hand  sides  of  Equation  (A4)  and  Equation  (3), 


2  Refer  to  the  Appendix. 
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we  have 


g"p(*»-g"r*.  (5) 

Since  g';  are  not  arbitrary,  we  conclude  that 

p*  *  -  r* 

r*i  ^  1M 

for  all  values  of  i,  j,  k. 

Recursive  Property  of  P* 

We  now  impose  the  following  requirement  on  Ps*:  If  the  coordinates  are 
such  that  their  Laplacians  vanish,  i.e.,  V  2xk  =  0,  then  the  control  functions 
P(*  vanish  individually  for  all  indices  i,  j,  k. 

The  importance  of  this  requirement  becomes  clear  by  considering  Equa¬ 
tion  (5)  in  which  T,*  are  the  Christoffel  symbols.  Thus,  when  the  Laplacians 
of  the  coordinates  vanish,  the  right-hand  side  of  Equation  (5)  vanishes  as  an 
inner  sum,  but  the  left-hand  side  vanishes  due  to  the  imposed  restriction  of 
“  individual  vanishing”  of  Pf-  for  such  coordinates.  Denoting  the  coordinates 
which  satisfy  the  Laplace  equations  by  x[0),  we  then  have 

Vm  “  r*(0>  *  o,  g  (0)^(0)  **  o*  (6) 

We  now  consider  those  coordinates  which  satisfy  the  Poisson  equations.  In 
these  coordinates  it  is  of  interest  to  know  the  relation  between  the  successive 
P ’s  under  successive  coordinate  transformations.  To  this  end,  we  consider  two 
successive  allowable  transformations  denoted  by  x(‘m_  n  and  x(‘m),  with  x‘0)  as 
those  coordinates  which  satisfy  the  Laplace  equation.  Thus,  in  the  r‘m_  n 
coordinates,  Equation  (4)  is 


a2r 


■"1  dxlrn-l)dx’m-l)  +  P,'<n 


dr 


'3*?-n 


=  0, 


(7) 


and  in  the  coordinates 


a2r 


dr 


dxln)dxU 


+  PP 


=  0. 


(8) 


(tn  )  , 
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Using  the  transformation  law  (AT)  in  the  form 


S(L- 1) 


dx?m)  3x"m)  g(m) 


(9) 


and  the  chain  rule  of  partiai  differentiation  of  r  in  Equation  (7),  and 
comparing  the  resulting  equation  with  Equation  (8),  we  get 


pp  = 


5*r. 


m) 


m) 


r  Arrt(  m  —  1) 


3xf 


dx*. 


m  —  1) 


0_n 

{m-  l)  °X(m-l) 


dXln» 


(10) 


establishing  a  relation  between  the  successive  P ’s.  This  is  how  Pk  transform 
from  one  coordinate  system  to  the  other.  Also,  using  Equation  (A9)  in 
Equation  (10),  we  have 


pp 


-IT U  +  te'-o- 


- 1)  l\n(im-  l)) 


3*fm)  3x(nm_u 


9*tm- 1)  ^xl 


m) 


(ID 


m) 


Equation  (11)  establishes  a  relation  between  the  P’s  and  F’s.  In  particular, 
for  m  =*  1,  Equation  (10)  reduces  to  Equation  (63)  of  Reference  [3]- 

From  the  recursive  relation  (10)  it  is  a  straightforward  matter  to  show  that 
the  transformation  x|0)  -*  x'w  x)2>  of  P*  is  the  same  as  the  transformation 
X(0)  *(V  The  use  of  the  chain  rule  for 


3x/i,  ’  ax/!,  5x*m 


and  of  the  formula 

aVo>  _  £%> 

ax/„  dx^  dx'[0)  dx‘0)  dx‘u  dx‘n)  9xfl} 


gives  the  required  result. 


Other  Poisson  Systems 

It  is  also  possible  to  have  equations  in  which  the  dependent  variables  are 
non-Cartesian,  e.g.,  cylindrical  or  spherical  coordinates.  In  a  paper  by  Ghia 
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et  al.  [13],  the  computational  coordinates  are  assumed  to  satisfy  a  Poisson 
,  system  in  the  cylindrical  coordinates. 

'  To  generalize  this  concept,  let  x*  be  a  coordinate  system  on  which  the 
coordinate  system  f  ’  (e.g.,  cylindrical)  is  to  be  generated.  From  Equation 

(A4), 


V2*’-  -grip*  (12)  •  ; 

and 

vV  -  -g<>r‘-g‘jp* 

=  P*  (say).  (13) 

In  Equation  (12)  both  the  metric  coefficients  and  the  Christoffel  symbols  are 
already  known  in  terms  of  the  x1  coordinate  system.  [If  x*  are  the  Cartesian 
coordinates,  then  Equation  (12)  is  an  identity.]  From  Equation  (A3)  we  then 
have 


7  *  i  - 


d*x’  u  dx* 

_  +.  p* - „  _  s'i r* 

i+  a**  g  d’ 


dx'dx 


(14) 


which  are  the  transformed  equations  in  the  computational  space.  By  using 
Equation  (A10)  we  can  also  write  Equation  (14)  as 


d2*’  „  .dx’  „  _ 


(15) 


where 


/  =  det 


Recently  Fujii  [14]  used  a  Cartesian-type  Poisson  system  between  two 
curvilinear  coordinate  systems.  Let  xJ  and  f ’  be  two  curvilinear  systems. 
Consider  the  Cartesian-type  Poisson  system 
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where  x'  is  the  curvilinear  system  on  which  the  x”  -system  (e.g.,  spherical)  is 
^to  be  generated.  To  obtain  the  inversion  of  Equation  (18),  we  take  the  inner 
product  of  Equation  (All)  with  dx'/dx'  to  obtain 

d *f‘  3x‘  dx>  3h'  dx‘ 

dx‘  dx>  IP  IP  =  ”  dxk  din  1? ’ 


Setting  k  =  n  l,  i  =»  n  =  2,  k  =■  n  =  3  and  adding,  we  get 


so  that 


where 


Cf' 


d2? 

dx*dx‘ 


£V 

5xr7  dx*dx’' 


C'l 


d2x' 
dx'  dx> 


(18) 


,  ^  dx’  dx> 

c’i=  r - 

7  di>  d£> 


Equation  (18)  is  the  inverted  form  of  Equation  (16).  Similarly  to  the  case  of 
actual  Poisson  system,  the  functions  Pr  are  again  the  arbitrary  specified 
functions.  Alternatively,  we  can  also  write 


P'-C'P,', 

and  then  are  the  arbitrary  control  functions. 

Beltramiuns  as  Grid  Generators  in  Curved  Surfaces 

The  problem  of  generating  curvilinear  coordinates  in  a  surface  has  im¬ 
portant  applications  in  many  branches  of  engineering.  In  this  paper  we  are 
interested  only  in  one  differential  model  which  generates  coordinates  in 
surfaces  embedded  in  fl  3  and  which  reduces  to  Equation  (2)  in  a  natural  wav 
when  the  surface  degenerates  into  a  plane.  To  achieve  this  aim,  we  consider 
the  formulae  of  Gauss  [Equation  (A12)J  for  a  surface  v  =*  const.  Inner 
multiplication  of  Equation  (A12)  by  C,gaS  yields 

Dr  +  G.(A'j'x*)r  4  =  n(''R,  (19a) 
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where 


D“G,gafidafi,  n(,> «  unit  surface  normal, 

C,  “  gaaSss  "  (s=^)2>  r,  a,  0  cyclic,  (19b) 

R~o,g^bafi~(k\'>  +  k\v)c„ 

K*  are  the  coefficients  of  the  second  fundamental  form,  and  k{’\  are  the 
principal  curvatures  at  a  point  of  the  surface  v  —  const.  The  Beltramians 
A1,’!3  are  defined  by  Equation  (A14).  It  is  readily  seen  that  when  the  surface 
v  —  const  degenerates  into  a  plane,  then  R  —  0  and  Az  becomes  the  Laplace 
operator. 

Equation  (19a)  can  be  used  either  for  the  generation  of  coordinates  in  a 
given  surface  or  for  the  generation  of  3D  spatial  coordinates  between  two 
given  surfaces.  In  the  latter  case  it  has  been  shown  [10]  that 

fl«C,g“*  X<'>«n<'>.r„  (19c) 

where  are  the  3-space  Christoffel  symbols  and  x'  is  the  transverse 
coordinate. 

Surface  Coordinates 

The  surface-oriented  generating  system  of  equations,  with  the  option  of 
arbitrary  coordinate  control,  is  now  obtained  by  putting  suitable  restrictions 
on  the  Beltramians  appearing  in  Equation  (19a).  Similar  to  Equation  (3),  a 
general  specification  of  the  Beltramians  is 

A^x’-g^V  (20) 

where  satisfy  all  the  properties  stated  earlier,  including  Equations  (10) 
and  (11)  with  T  replaced  by  T. 

To  be  specific,  we  take  the  surface  x3  =  const  as  the  given  surface  and 
x1  —  £,  x2  =  tj.  Then  Equation  (19a)  becomes 

Lr  =  n,3’fi ,  (21a) 


where 

~  +  +  F<?(  +  Q^rr 

^  =  gijf’u  ~  2gl2/,1j guFn),,  (21b) 

Q  “  S-2S^U  ~  |2 F, 2  +  g(1Fj|. 
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3.  SURFACE-GRID  EQUATIONS  FROM  THE  LAPLACIANS 


In  this  section  an  attempt  is  made  to  establish  a  relation  between  the 
surface  equations  derived  by  using  the  formulae  of  Gauss  [Equation  (19a)] 
and  the  surface  equations  as  derived  from  the  inversions  of  the  3D  Laplacians. 

Let  £,  7),  £  be  a  general  curvilinear  coordinate  system  in  R3  such  that  £,  tj 
form  the  coordinates  in  the  surface  f » const,  with  f  as  the  transverse 
coordinate.  Starting  from  Equation  (B.3)  of  Reference  [1]  and  using  Equation 
(A2),  we  get  (refer  to  [11]  for  more  details) 


Qr 


A 


where  the  differential  operator  Q  is 


Q 


D  +  T“ 


dx “ 


-D  +  (S“  +  C3  A4z“) 


(22a) 


(22b) 

(22c) 


In  (22b,  c)  the  differential  operator  D  defined  in  (19b)  is 

D”g»d{|-2g122{„  +  g  u3„.  (22d) 

and 

T"  ”  2g [jrj^  —  gj^rj)  —  gur^,  (22e) 

sa  =  2g  12(  r“  -  T  ji )  -  g*(  rr,  -  Tf, )  -  g „(  r&  -  T£ ).  (22f ) 


where  a  =  1,2.  Further 


A  =  n-rf,  (22g) 

n  being  the  unit  surface  normal  vector  on  f  =  const. 

It  must  be  emphasized  that  T'k  and  are  the  3-space  and  2-space 
Christoffel  symbols  of  the  second  idnd  respectively,  and  in  general  they  are 
not  equal  to  each  other  at  f  =  const.  Thus,  a  comparison  of  Equations  ( 19a) 
and  (22a)  shows  that  the  two  sets  of  surface  generators  are  entirely  different. 
However,  it  has  been  shown  below  that  Equations  (19a)  and  (22a)  become 
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:  r.?  v  v*;:  ^ 

exactly  the  same  equations  under  the  following  conditions: 

v  (i)  When  the  surface  degenerates  into  a  plane,  in  which  case  £,  t  are  the 

curvilinear  coordinates  (of  any  nature)  in  a  plane. 

(ii)  When  £,  tj  are  any  general  coordinates  in  a  surface,  but  $  (the 
transverse  coordinate)  is  orthogonal  to  the  surface. 

Case  (i)  is  patently  straightforward.  For  case  (ii),  noting  that  f  is  orthogo¬ 
nal  to  the  surface,  we  have 

813  ”  823  “ 

(23a) 

-• 

and 

* 'r  ‘i 

rr  ,  1 — 

n  ”  ^  *  yg33  • 

(23b) 

Under  the  conditions  (23a)  it  is  easy  to  show  that 

r;r-T/y, 

where  all  Greek  indices  assume  only  two  values  (here  1,2).  Thus  Sa  —  0,  and 
using  (23b)  we  find  that  Equations  (19a)  and  (23a)  become  the  same 
equations  under  the  condition  of  orthogonality  of  the  f-coordinate  to  the 
surface  f  —  const.  Further,  under  this  condition  the  Beltramians  and  the 
Laplacians  of  £  and  if  are  related  through  the  following  equations: 

V2f-A2£-  — , 

S33 

(24a) 

2  , 

A,7) - . 

833 

(24b) 

Further, 

v-f  =  , —  , 

833  V  833 

(24c) 

k\3)  +  fc'f  ■ 


_J _ 3C, 

2^G3  dS  ’ 


(24d) 


where 
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and  (referring  to  (19b)J 


C3-gn«a-(«i a)*- 

4  SURFACE  GENERATING  EQUATIONS  BASED  ON  A 
VARIATIONAL  PRINCIPLE 

The  use  of  a  variational  principle  in  the  generation  of  both  2D  and  3D 
grids  has  already  been  considered  by  Brackbill  and  Saltzman  [15,  16]  and 
Thompson  and  Warsi  [17].  In  this  section  we  shall  consider  only  the  surface 
grid  generation  problem  based  on  the  use  of  a  variational  principle.  In  essence 
the  following  analysis  is  a  unified  approach  to  both  the  plane  2D  and  surface 
2D  cases  from  a  variational  viewpoint. 

Let  r1,  x2  be  the  coordinates  in  a  surface.  Consider  the  surface  functional 

/- fjcl+dx'dx2,  (25) 

where  C3~  gngn-  (&n)2  and  ^  is  a  specified  function.  The  condition 
SI  -  0  then  leads  one  to  the  Euler-Lagrange  equations  (using  the  summation 
convention) 


3  I  Jet)  -  iSiii-o 


.0 


(26) 


where  xr  (r  -  1,2,3)  are  the  rectangular  Cartesian  coordinates,  xB  are  the 
curvilinear  surface  coordinates,  and 


dx,  dh, 

X,B~J?'  x,a*~  IFl?' 

From  Equation  (26),  it  is  a  purely  algebraic  problem  to  show  that 

y'c;  d  l  1  dC1 


where 
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Let  now  be  a  function  of  x,fl.  Then  expanding  Equation  (26)  and  using 
Equation  (27),  we  get 


Lx,  =  -  M, 


where 


(28a) 


dC3  d<f>  d<t>  dC3 

c,J_ 

d  <f>  \ 

dxr  J  dx 3  dxr g  dx“  l 

<p  dx& 

,dx'.0l 

<f> 


On  the  other  hand,  if  £  is  taken  as 

*  =  J7C3, 

where  F  is  still  a  function  of  x,  then  in  place  of  (28b),  we  have 


1  (  dC3  dF  dC3  dF 
M  “  —  I  — — —  ;  + 


2E  \  dxr  „  dx 0  dx*  dxr 


i  dc3  dc3  c3  a 

lh*~T  IT* 


dF 


dx 


(28b) 


(29a) 


(29b) 


The  generating  system  (28a)  with  L  defined  in  (27b)  is  similar  to  Equation 
( 19a).  However,  the  selection  of  the  form  of  the  function  <f>  or  F  which  yields 
the  right-hand  side  of  Equation  (19a)  seems  to  be  a  difficult  task.  One  simple 
case  is  when  <(>  —  1.  In  this  case  the  minimization  of  /  implies  [from  (28a)] 

Lx,  **  0, 


and  these  are  the  equations  for  a  minimal  surface.  The  form  (29a)  is  of 
interest  because  the  choice 


F  S 1 1  +  6  22  - 


(30a) 


or 


<P  ”  g 


n 


(30b) 


is  equivalent  to  the  "smoothness'’  problem  in  2D  plane  coordinates,  as  has 
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been  shown  in  Reference  [I5J.  It  must,  however,  be  stated  that  “smoothness” 
of  coordinates  in  a  2D  plane  problem  is  due  to  the  satisfaction  of  the  Laplace 
equation.  No  such  criterion  is  obvious  on  using  either  of  the  equations  (30)  in 
(28a). 


5.  ANALYTICAL  COMPARISONS  AND  CONCLUSIONS 

Based  on  the  foregoing  analysis,  we  conclude  as  follows: 

(a)  The  Laplace-Poisson  system  for  2D  regions  is  optimal,  since  its  inver¬ 
sion  coincides  exactly  with  the  Gaussian  equations  in  a  plane.  (Note:  In  a 
plane  the  right-hand  side  of  Equation  (19a)  is  zero  and  the  Beltramians 
become  the  Laplacians.] 

(b)  The  inversion  of  the  Laplace-Poisson  system  for  3D  regions  for  a 
constant  coordinate  value,  viz.,  for  a  surface,  does  not  reduce  to  the  Gaussian 
equations  (cl.  Section  3)  except  when  the  transverse  coordinate  is  taken  as 
orthogonal  to  the  surface  [cf.  Equation  (22a)].  This  implies  that  the  Laplace- 
Poisson  system  in  3D  regions  with  three  nonorthogonal  coordinates  is  not 
optimal,  though  it  is  a  valid  system.  The  extra  terms  (22f)  should  somehow  be 
managed,  and  in  practice,  they  are  taken  as  part  of  the  arbitrary  specified 
control  functions.  This  means  that  the  generated  coordinates  will  assume  a 
distribution  which  may  not  be  to  one's  desire.  It  must,  however,  be  again 
emphasized  that  the  terms  (22f)  vanish  when  the  transverse  coordinate  is 
orthogonal  to  the  surface.  (In  Reference  [18],  the  author  had  to  make  other 
assumptions  besides  orthogonality.) 

(c)  For  the  elliptic  system  described  by  Equation  (15)  the  conclusions 
discussed  in  (b)  hold  good. 

(d)  The  generating  system  described  by  Equation  (16)  is  a  Cartesian-type 
Poisson  system.  It  looks  difficult  to  assess  its  optimality  in  relation  to  the 
Caussian  equations. 

(e)  The  equations  derived  from  the  variational  principle,  viz.,  (28),  are 
nearly  optimal,  though  it  looks  difficult  to  find  the  appropriate  function  <?> 
which  makes  the  right  hand  side  of  (28a)  the  same  as  that  of  Equation  (19a). 


APPENDIX 

The  following  formulae  have  been  used  in  the  main  text  of  the  paper  and 
can  be  found  in  any  standard  text  on  tensors,  e.g.,  [19],  [20], 

In  an  Euclidean  space  R”  (though  here  we  are  concerned  only  with  R1  or 
R1),  endowed  with  the  curvilinear  coordinates  x‘,  the  first  partial  denvatives 
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of  the  covariant  basis  vectors  a,  are  expressible  as  linear  functions  of  a,: 


da, 

dx' 


r*j 


(Al) 


where  F,k;  are  the  Christoffel  symbols  of  the  second  kind.  Here,  and  in  what 
follows,  a  repeated  lower  and  upper  index  will  always  imply  summation  over 
the  range  of  index  values.  Equation  (Al)  can  alternatively  be  written  as 

r./=  T,kj r.k  •  (A2) 


where  a  comma  preceding  an  index  implies  a  partial  derivative,  and  r  is  the 
position  vector,  viz., 

r  =  (x,y,z)or(x1,x2,x3). 

The  Laplacian  of  a  scalar  <j>  is  given  by 

(A3) 

From  (A3),  we  have 

V2x*»-g,'r*  (A4) 

The  quantities  g ’’  and  g,j  are  respectively  the  contravariant  and  covariant 
components  of  the  metric  tensor,  and  the  two  are  related  as 


gpngr*  gpfg  nr 


where  (i.p.r)  and  ( j,n,s )  are  to  be  taken  in  the  cyclic  permutations  of 
(1,2,3),  and 


g“det(gt>). 

Let  x‘  and  x'  be  two  allowable  coordinate  systems  in  a  Euclidean  space 
such  that  each  of  the  functions  x'  =  <t>(x')  and  x’  =  <//'(x;)  define  a  transfor¬ 
mation  with 


/  =*  det 


*  0, 


/  =  det 


*0. 


H 
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On  transformation  of  coordinates  from  x'  to  x‘,  we  have 

_  8  £ 

S  “  ji<  8  "  jz  > 

dx'  dx’ 
gp"  =  1? 

dxp  dxn 

» - g'i 

s  3*'  3xj B  ’ 

_n  ^  dx”  3x‘  dx>  3  V  dxp 

rf«  “  r«f  aP  +  aFdF  3x>  ’ 

d%p  int  3xp  _o3xr  3x' 

3xk3x"  “  ~  T^17klP' 


(A5) 

(A6) 

(A7) 

(A8) 

(A9) 


The  first  partial  derivatives  of  x‘  with  respect  to  x1  are  given  by 


dx1  “  ] 

where 


dxT  dxk  3xr  3xk 
C'  "  3?  V"  ~  fo"  dT'  ’ 


(A10) 


and  (i,  r,  n)  and  (;,  r,  k )  are  cyclic  permutations  of  (1,2,3). 

The  second  partial  derivatives  of  one  set  of  coordinates  are  related  with 
those  of  the  other  set  as 


3V  3V  dxr  9x'  dx> 

3r‘  3£"  dx'  dx’  IT’  JTk  3f"  ' 


(All) 


We  now  consider  a  surface  embedded  in  RJ.  All  tensor  indices  associated 
with  a  surface  will  be  denoted  by  Greek  letters  (except  the  letter  v ).  In 
contrast  to  Equation  (A2),  the  formula  of  Gauss  is 

r.s+nMbof  (A12) 
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where  are  the  surface  Christoffel  symbols  of  the  second  land,  baB  are  the 
coefficients  of  the  second  fundamental  form,  and  n(r)  is  the  unit  surface 
normal  on  the  surface  v  =  const.  Tlie  values  of  v  and  other  Greek  indices 
follow  the  following  scheme: 

v  =  1:  Greek  indices  a,  0,  etc.  assume  values  2  and  3. 

v  =  2:  Creek  indices  a,  0,  etc.  assume  values  3  and  1. 

v  =  3:  Greek  indices  a,  0,  etc.  assume  values  1  and  2. 

The  Beltramian  of  a  scalar  <f>  is  given  by 

-T „>.,)•  (A13) 

From  (A13),  we  have 

A(2''x4= -g^T„V  (A14) 

For  a  surface  the  formulae  (A5)-(A9)  are  equally  applicable  with  proper 
cho.ce  of  indices  and  replacing  T  by  T. 
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